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Abstract 

From  any  given  Iterated  Function  System,  a  small  set  of  balls  that  cover  the 
fractal  attractor  can  bo  simply  determined.  This  gives  a  priori  bounds  on 
the  region  of  space  in  which  the  attractor  may  be  constructed. 


As  Barnsley.  Demko  and  others  have  shown  [1.  2.  C.  4],  one  effective 
method  for  producing  fractal  shapes  (in  any  number  of  dimensions)  is  with 
Iterated  Function  Systems  (IFSs).  using  the  "Chaos  Game"  algorithm  (or 
some  deterministic  algorithm).  This  approach  has  been  used  for  producing 
naturalistic  shapes  [4],  finding  fractal  interpolants  to  given  data  {6.  p.  274] 
and  fractal  approximations  of  given  functions  [8].  and  even  for  visualizing 
arbitrary  discrete  sequences  [7|.  Indeed,  any  contractive  IFS  will  give  an 
attractor  (usually  of  fractal  dimension);  thus  it  is  possible  to  generate  IFSs  at 
random  to  explore  the  graphical  possibilities,  as  is  done  in  some  educational 
software  [5],  Similarly,  because  the  attractor  depends  continuously  on  the 
parameters  in  the  IFS  [l],  a  small  data  sets  from  any  source  could  be  encoded 
as  IFSs  for  visualization. 

In  implementing  the  IFS  method,  one  important  question  is  the  prediction 
u  priori  of  the  region  of  space  containing  the  fractal  attractor.  Without  such 
a  prediction,  one  could  only  approximately  estimate  the  spatial  extent  based 
on  calculating  several  points  of  the  attractor,  with  no  guarantee  that  these 
points  are  near  the  bounds.  If  as  a  result  the  portion  of  space  represented 
in  the  computation  of  the  attractor  is  too  small,  the  result  will  not  yield 
the  whole  attractor.  If  the  portion  of  space  represented  is  overly  large,  then 
much  computational  space  is  wasted,  reducing  the  effective  resolution  of  the 
computed  attractor.  Another  concern  is  when  the  space  has  natural  limits, 
for  example  the  space  of  colors  in  Red-Green-Blue-space  representable  on  a 
video  monitor  is  limited  (more  or  less)  to  a  unit  cube.  Then  the  question 
becomes  whether  the  attractor  (or  its  projection  onto  the  limited  dimensions) 
will  fit  in  the  space. 

Here  a  simple  algorithm  is  given  to  compute,  directly  from  the  IFS.  a  set 
of  balls  whose  union  contains  the  attractor  as  a  subset  (see  Figure  1).  The 
radii  of  the  balls  are  minimal  in  a  certain  restricted  sense.  This  gives  reliable 
bounds  on  the  region  of  space  that  must  be  considered  in  constructing  the 
fractal.  The  method  is  general,  independent  of  the  particular  space  and 
metric.  We  first  describe  the  set  of  balls,  then  show  how  to  compute  their 
radii  and  prove  that  the  algorithm  works,  and  lastly  give  a  detailed  example. 

An  IFS  consists  of  a  set  of  n  contraction  mappings  w,  :  X  -+  X  on 
a  metric  space  X  with  metric  cl  :  X  x  X  -»  TL.  (For  the  "Chaos  Game’’ 
algorithm,  probabilities  p,  are  associated  with  each  mapping;  this  idea  has 
been  extended  to  conditional  probabilities  [3].  Here  only  the  former  case  is 
considered,  where  the  attractor  is  independent  of  the  non-zero  p,.)  Assume 


1 


♦ 


that  for  each  contraction  w,  the  contractivity  ratio  0  <  s,  <  1  and  the  fixed 
point  x,  are  known,  where  by  definition  s,  satisfies  d(u',(x),  «•,( y))  <  s,d(x.y) 
for  all  points  x.y  €  X  and  x,  satisfies  x,  =  w,{x,).  The  action  of  the  IFS  H' 
on  a  set  S  of  points  in  X  is  defined  as 

H*(S)  =  U  Vi(S)  (lj 

.=i 

where  each  contraction  ir,  is  applied  to  the  set  5  in  a  pointwise  sense.  The 
attractor  A  is  the  set  of  points  in  X  satisfying 

A  =  W(A)  .  (2) 

That  is.  the  attractor  consists  of  n  smaller  "copies'"  of  itself. 

We  seek  to  cover  each  of  the  »  "copies"  with  a  closed  ball  B,  centered 
on  the  corresponding  fixed  point  x,.  so  the  radius  r,  must  be  chosen  large 
enough  that  B,  D  ie,(A).  Call  the  union  E  of  these  balls  the  "envelope."  in 
that  E  D  A  by  (2).  Then  relative  to  each  x,.  every  point  in  the  envelope  will 
be  within  a  distance  R,  =  niaXjfr/.j  +  ij).  where  d,}  =  d(x,.Xj).  because  for 
any  point  x  in  By  d(x,.x)  <  <l,j  +d(x;.x)  <  (ltJ  +  ry  Applying  w,  to  such  a 
point  x  will  give  an  image  point  y.  where  <l(x,.y)  <  $,d{x,.x)  <  Hence 
if  the  radii  r,  are  chosen  to  satisfy 


r,  =  s,  maxfc/,,  +  /’, )  (3) 

for  i.j  =  1  . . .  v  then  B,  will  contain  the  image  ic, ( JET )  of  the  envelope  and  so 
E  will  contain  its  own  image  under  the  IFS: 

EDW(E)  .  (4) 


Iterating  the  IFS  from  any  starting  set  (E  in  particular)  yields  a  sequence  of 
sets  that  converges  to  the  attractor.  Since  (4)  implies  E  D  U'a'(jF)  for  any 
positive  integer  A\  the  envelope  E.  subject  to  (3).  does  indeed  contain  the 
attractor.  But  how  can  the  r,  be  calculated  from  (3)? 

When  n  =  2  the  radii  can  be  determined  algebraically.  Solving  the  pair 
of  equations  (3)  gives: 


r  i 


r-2 


51 (1  A  s2) 

1  -  SiS-2 

5a(l  +  s\) 

1  —  S  i  S) 


d\ 2 


dV2 


(5) 
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* 
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But  for  n  >  2  there  is  apparently  no  closed-form  general  solution,  and  the  r, 
must  be  found  algorithmically.  (If  one  wished  to  minimize  calculating  at  the 
expense  of  overestimating,  one  could  use  r,  =  dmaTsmar/(  1  —  i.e..  the 

envelope  as  if  all  fixed  points  were  equidistant  at  the  maximum  separation 
and  with  all  s,  equal  to  the  largest.) 

A  natural  approach  for  n  >  2  is  to  start  with  the  pairwise  estimates 


,()i 


1  -  S,Sj 


d,j 


max  r,. 

j*' 


(C) 


but  in  most  cases  the  r,(1)  will  not  satisfy  (3).  The  exceptional  case  is  when 
r,j  =  for  every  i.j  ^  i.k  #  i.  i.e..  when  for  each  ball  all  the  pairwise 
estimates  for  that  ball  give  the  same  size.  (This  case  is  not  always  apparent 
from  the  attractor:  Figure  1  shows  such  an  example.)  Otherwise,  some  of 
the  rj11  will  be  too  small  to  contain  some  images  )  of  the  other  balls. 

Then  the  obvious  iterative  scheme  to  try  is 

r,U+1)  =  max.s,b/,,  +  rlk> )  .  /.j  =  l...n  .  (S) 

Because  this  approach  never  overestimates  the  radii  (r[k>  <  r,)  and  the  iter¬ 
ates  are  nondecreasing.  the  algorithm  must  converge.  What  is  not  so  obvious 
is  that  this  process  always  sit  coeds  in  at  most  n  —  1  iterations,  as  shown  be¬ 
low.  (This  would  not  be  true  without  (7).)  In  fact,  there  is  a  direct  algorithm 
(not  iterative)  that  is  more  efficient. 

The  key  idea  that  allows  the  direct  algorithm  is  that  the  distances  d,j  can 
be  rescaled  to  account  for  the  eontractivities  s,.  and  the  scaled  distances  DtJ 
can  be  used  to  order  the  contractions  tc,.  Let 


__  (1  +  s,)(l  +  Sj)  , 

I J  —  •,  ^1> 

1  -  S,Sj 


J,  j  =  1 . .  .n 


(9) 


(While  D,j  is  clearly  symmetric  and  non-negative,  it  is  not  a  metric  because 
it  doesn't  satisfy  the  triangle  inequality.)  Now  reorder  (and  relabel)  the  w, 
bv  decreasing  maximum  scaled  distance,  so  that 


i  <  j  =>  max D,l  >  max Dji  .  i.j.k.I  =  !...»?  (10) 


3 


In  the  new  order,  use  the  pairwise  formula  (5)  for  the  first  two  radii.  Then 
proceed  in  order  based  on  the  previous  results,  letting 


r,  =  s,  max((/„  -f  r . )  t-  3.. 


(lit 


This  is  the  direct  algorithm,  which,  as  shown  below,  solves  (3):  an  imple¬ 
mentation  in  the  C  programming  language  is  given  in  the  Appendix. 

First  note  that  in  the  exceptional  case  mentioned  above  (after  7).  all  the 
D,j  (for  i  ^  j)  are  equal.  In  this  case,  the  direct  algorithm  will  obtain  the 
correct  r,  regardless  of  the  order  in  which  they  are  computed.  If  not  all  the 
D,j  are  equal,  then  some  of  the  r,  will  need  to  be  larger  than  the  pairwix 
estimates  r 


and  hence  larger  than  the  r 


r,  > 


s, 


1  +  A, 


tJ  in  (6).  so  in  general  (3)  implies 

D,j  (12) 


for  i.j  =  l...n. 

For  the  general  case,  the  proof  is  by  induction,  showing  that  each  new  i , 
computed  lequires  no  adjustment  of  those  previously  computed.  Clearly  /  j. 
70  from  (5)  satisfy  (3)  for  the  subset  i.j  =  1,2.  Now  assume  the  first  m  -  1 
radii,  in  the  order  (10).  satisfy  (3).  and  hence  (12).  for  i.j  —  1  . . .  w  —  1 . 
Choose  r„,  by  (11).  and  let  k  be  the  value  of  the  index  j  in  (11)  for  which 
the  maximum  is  achieved.  Then  by  (12)  and  the  ordering  (10) 
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Sh 


1  +  At 

Algebraic  manipulation  of  (13)  gives 


Di:, 


^k  ^2  A  #,[ )  1  4“  A,„  4“  A , , ,  /  r  —  A  ( (i  kru  4"  1  rn  ) 


(13) 


1-4 ) 


so  the  new  /■„,  requires  no  alteration  of  rk. 
Similarly,  for  /  ^  k.  i  <  in 

s 


r,  > 


i 

1  4-  s, 

Si  ( 1  4-  sni ) 

1  -  A ... 


dtr 


(15) 


since  s,  <  1.  Combining  (15)  with 


r,  >  s,(d,k  4-  rk) 


TG) 


« 
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yields 


^  ^i[(l  d"  Srn)(l,m  +  Sm(d,k  •+■  jq.)] 

^  5,[d,m  "f*  Sm  (dkm  4“  7q.  )j  ”  $,  {dtni  “4“  rm  )  ill) 

by  the  triangle  inequality.  Hence  the  new  rm  requires  no  adjustment  of  any  r, 
for  i  <  m.  and  (3)  is  satisfied  for  i,j  =  1  .  . .  m:  this  completes  the  proof.  Note 
also  that  because  of  how  the  direct  algorithm  works,  the  iterative  algorithm 
will  compute  at  least  two  of  the  r,  (r, .  r2  in  the  order  (10))  in  the  initial  step, 
and  will  find  at  least  one  of  the  other  r,  at  each  successive  step,  and  so  can 
take  at  most  —  1  iterations  to  arrive  at  the  answer. 

Can  radii  smaller  than  these  r,  be  used  and  still  have  the  D,  covet  the 
attractor?  For  any  particular  IFS.  the  answer  is  probably  yes  (as  illustrated 
in  Figure  2).  The  approach  given  here  uses  only  minimal  information  about 
the  IFS:  the  ordering  of  the  it',  based  on  the  maximal  scaled  distan  •<*>:  tin- 

and  for  each  i  one  determining  distance  d,k  (where  the  maximum  in  ill 
is  achieved).  Using  more  infoimation  it  may  be  possible  to  reduce  the  si/e 
of  the  D,.  But  if  one  considers  the  set  of  all  the  IFSs  for  which  the  diteet 
algorithm  yields  the  same  r,  in  the  same  way  (i.e..  same  ordered  *,  and  same 
>i  -  1  determining  distances),  then  tin*  r,  are  minimal  for  that  set  of  attractors 
(see  Figure  3).  In  fact,  one  can  construct  one  member  of  that  set  such  that 
each  image  w,[A)  of  the  attractor  includes  a  point  at  a  distance  r,  from  the 
fixed  point  r,. 

To  construct  this  IFS.  let  X  =  TZ  with  the  Euclidean  metric  r/(.i  . .//)  = 
|.r  -  y\.  Let  te,{.v)  =  r  +  (1  +  s,).r,.  and  let  .iq  =  0.  ,;q  =  d say.  Then 
the  attractor  includes  the  extremal  points  ,j  el  =  jq  —  /q.  .if  ,  =  .i  j  +  since 
.r,i  =  uq  (.r(  > )  and  .(  ,•>  =  imbUib  Place  each  succeding  .r„,  at  the  determining 
distance  d„lk  from  the  determining  point  ,jq.  in  the  opposite  direction  from 
.rck.  (Figure  4  illustrates  the  construction.)  Then  the  attractor  will  include 
•i 'em  =  ±  >■„,  =  ie„,  (./>*.).  Thus  for  this  one-dimensional  attractor  A.  each 

image  u',(A)  will  include  a  point  (,rtl)  a  distance  r,  from  r,.  so  no  smaller  r, 
would  suffice. 

As  an  example,  consider  the  now  familiar  black  spleenwort  fern  fractal  of 
[4].  The  IFS  for  the  fern  (in  two  dimensions)  consists  of  affine  contractions, 
each  of  which  has  the  form 


o 


or  rap;.'  ompactly 


U',(x)  =  A/,x  +  b, 


(19) 


where  A/,  is  the  matrix  and  b,  is  the  offset  vector.  The  various  constants  are 
given  in  [-1],  but  in  terms  of  scaling  and  rotating  each  axis,  using  p.  q.  B.  o. 
where  a  —  pc osB.  b  =  — c/sino,  c  =  psin#.  d  —  qcoso.  The  following  table 


is  adapted  from  4, 

p.  197  /  : 

Map 

Scalings 

Rotations 

Translations 

i 

P. 

<1> 

o, 

o, 

e, 

/. 

1 

o 

0.10 

0 

0 

0 

0 

2 

0.83 

0.83 

-2.5 

-2  5 

0 

l.G 

3 

0.3 

0.34 

49 

49 

0 

1.0 

4 

0.3 

0.37 

120 

-50 

0 

0.44 

where  angles  are  given  in  degrees.  The  contractivity  ratio  s.  for  an  affine 
map  u1,  is  the  largest  singular  value  of  the  matrix  A/,.  In  the  first  three  maps 
above,  both  axes  rotate  together,  and  so  s,  is  the  larger  of  p,.  q, .  In  a.,,  the 
differential  rotation  causes  a  skewing  effect,  and  the  singular  values  of  A/< 
must  be  found.  The  simplest  way  for  a  real  2x2  matrix  is  first  factoi  out 
a  pure  rotation  to  give  a  symmetric  matrix  (5).  then  diagonal!”  it  to  find 
its  eigenvectors  (Aj.  A_>)  as  shown  below: 


o 


3 


At  0  \ 

0  A,  ) 


arctan 


c-h 
a  +  (I 


(coso  sino 
—  sin  n  cos  n 


1 

-  arctan 
2 


(  cos  i  sin  i 
y  —  sin  3  cos  3 


A/ 


S 


cos  3 
sin  3 


—  sin  3  A 
cos  3  J 


(20 
(21) 
1 22 . 
(23/ 


Then  ,s  =  max(|Aj  |.  |A>|).  This  approach  also  has  a  nice  geometrical  inter¬ 
pretation:  the  effect  of  multiplying  a  vector  x  by  A/  is  to  take  components 
of  x  in  the  eigenvector  directions,  which  are  at  an  angle  3  relative  to  the 
coordinate  axes,  scale  each  component  by  the  corresponding  A.  and  rotate 
the  resulting  vector  by  o. 

Proceeding  as  above  gives  s4  =  0.379.  The  fixed  point  x,  for  each  map 


0 


can  be  found  bv  solving 


(I  —  M,)  x,  =  b,  (24  5 

This  gives  the  necessary  starting  information,  summarized  as  ’  Input"  in  the 
table  below.  For  affine  maps  in  higher  dimensions,  the  contractivity  ratios 
are  found  by  singular  value  decomposition,  but  for  nonlinear  maps  the  ratios 
and  fixed  points  may  be  more  difficult  to  find. 

Running  the  direct  algorithm  program  ("envelope")  from  the  Appendix 
on  this  input  gives  the  following  results: 

envelope  4 

Map  1.  Enter  scale,  x,  y:  .16  0  0 
Map  1:  s  =  0.160000,  x  =  0.000000,  y  =  0.000000 
Map  2.  Enter  scale,  x,  y:  .88  2.45967  10.004734 
Map  2:  s  =  0.850000,  x  =  2.459670,  y  =  10.004734 
Map  3.  Enter  scale,  x,  y:  .34  -0.601889  1.883961 
Map  3:  s  =  0.340000,  x  =  -0.601889,  y  =  1.883961 
Map  4.  Enter  scale,  x,  y:  0.379216  0.155336  0.630251 
Map  4:  s  =  0.379216,  x  =  0.155336,  y  =  0.630251 
radii  in  sorted  order[orig  order] (sorted  link): 
rl[2](->2):  16.700212 
r2 [4] (->1 ) :  9.993765 
r3 [3] (->1 ) :  8.628835 
r 4 [1] (->1) :  4.320459 

7. 


These  results  are  illustrated  in  Figure  5.  and  tabulated  below  (including  the 
D  values  used  in  re-ordering  and  the  determining  distances): 


Input 

Output 

i 

a, 

■f, 

.  .. 

D-order 

D 

max 

(I 

d<  t 

1 

0.16 

0 

0 

4 

Dn. 

25.59 

dl2: 

10.30 

4.32 

2 

085 

2.460 

10.005 

1 

D2A: 

36.35 

d 24 •' 

9.65 

16.70 

3 

0.34 

-0.602 

1.884 

3 

D22'- 

30.26 

^32' 

8.68 

8.G3 

4 

0.379 

0.155 

0.630 

2 

Da2: 

36.35 

dA2' 

9.65 

9.99 

(While  in  this  example  the  determining  distance  for  each  map  derives  from 
the  same  pair  that  gives  the  maximum  D ,  that  is  not  always  the  case.)  Then. 
If  we  had  no  idea  how  big  the  fern  attractor  was.  we  could  use  a  computational 
space  extending  from  .r„„„  =  .r2  -  /•>  =  —14.24  to  xmaT  =  ,i  >  +  r->  ~  19. 1C 


and  ymtn  —  yA  —  r4  =  — 9.3G  to  ymaz  =  V2  +  r2  —  26.70  to  contain  the  entire 
envelope.  As  it  turns  out.  this  is  far  more  space  than  necessary  for  the  fern 
itself,  but  there  are  many  other  IFSs.  equivalent  as  far  as  the  direct  algorithm 
is  concerned,  with  much  larger  attractors  (e  g.,  what  if  Q2  =  o>  =  177.5 
instead). 

To  summarize,  given  any  IFS  (along  with  the  contractivities  and  fixed 
points  of  each  of  its  constituent  contraction  mappings),  an  envelope  can  be 
constructed  of  one  ball  for  each  map.  centered  on  the  corresponding  fixed 
point.  (In  the  case  of  affine  maps  in  two  dimensions,  an  explicit  proceduo 
for  finding  the  contractivities  and  fixed  points  was  given.)  We  have  proven 
that  the  radii  of  the  balls  can  be  calculated  by  a  simple  algorithm  (direct 
or  iterative)  such  that  the  envelope  covers  the  entire  attractor  .  The  spatial 
extent  of  the  envelope  thus  gives  a  reliable  bound  on  that  of  the  attractor.  (In 
addition,  if  the  balls  are  disjoint,  the  attractor  is  totally  disconnected. )  While 
the  radii  found  by  the  direct  algorithm  may  not  be  minimal  for  the  particular 
IFS.  they  are  minimal  for  the  set  of  all  IFSs  with  oouivalent  information  m 
the  sense  described  above). 
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APPENDIX 

f* 

The  following  program  in  C  implements  the  direct  algorithm 
for  determining  the  "envelope"  of  an  attractor  of  an  Iterated 
Function  System  on  R"2,  given  the  contractivities  si  and  the 
fixed  points  (xi.yi). 

written  by  David  Canright,  March  1993. 

*/ 

•include  <stdio.h> 

•include  <math.h> 

int  npts,  i,  j,  m,  n,  index [64],  1 ink [64]; 
double  d[64][64],  x[64],  y[64],  s[64],  r[64],  Dmax[64]  , 
t,  tmax,  dx,  dy; 

int  input (int  arge,  char  *argv[]>; 

mainCint  arge,  char  *argv[])  { 

npts  =  input(argc,  argv);  /*  get  si,  xi,  yi  */ 

/ *  compute  distances  dij  b  maximal  scaled  Dij  */ 
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for  (i  =  1;  i  <=  npts;  i++)  Dmaxti]  =0.,  » 

for  Ci  =  1;  i  <=  npts;  i++)  { 

for  (j  -  i+1;  j  <=  npts;  j++)  { 

dx  =  x[i]-x[j];  dy  =  y[i]-y[j];  « 

d[i]  [j]  =  d[j] [i]  =  t  =  sqrt (dx*dx+dy*dy) ; 
t  =  (l.+s[i])*(l.+s[j])/(l.-s[i]*s[j3)  *  t; 
if  (t  >  Dmaxfi])  Dmax[a]  =  t; 

if  (t  >  Dmaxfj])  Dmax[j]  =  t; 

} 

> 

/*  Sort  by  scaled  distances;  index  points  to  old  order  */ 
index  [1]  =  1; 

for  (i  =  2;  l  <=  npts;  i++)  { 

for  (m  =  i;  m  >  1  &&  Dmaxfi]  >  Dmax [index [m-1] ] ;  m — ) 
indexfm]  =  index[m-l]; 
index [m]  =  i ; 

> 

/*  Direct  algorithm;  link  points  to  determining  distance  */ 
i  =  index[l];  j  =  index [2] ;  linkfl]  =  2;  link[2]  =  1; 
r[l]  =  (  s[i]/(l  ,+s[i] )  )  *  Dmaxfi]; 
r  [2]  =  (  s[j]/(l  .+s[j] )  )  *  Dmax  [j]  ; 
for  (m=3;  m  <=  npts:  m++)  { 
i  =  index [m] ; 
tmax  =  0 . ; 

for  (n  =  1 ;  n  <  m;  n++)  { 
j  =  index [n] ; 

if  C  (t  =  dfi]  [j]  +  rfn])  >  tmax  ) 

{  tmax  =  t;  linkfm]  =  n;  > 

> 

rfm]  =  s Ci]  *  tmax; 

> 

printf ("radii  in  sorted  orderforig  order] (sorted  link):\n"); 
for  (i  *  1;  i  <=  npts;  i++) 

printf  ("  r’/.dfXd]  (->’/,d) :  ’/.f\n"  ,  i ,  indexfi]  .linkfi]  ,r[i]  ); 
return(O);  • 
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/*  Input  function:  gets  si,  xi,  yi,  or  uses  random  numbers  */ 

/*  optional  arguments:  number  of  maps,  seed  for  random  */ 

*  int  input (int  argc,  char  *argv[])  { 

double  norm; 

char  line  [81],  getmore  =  1; 

npts  =3;  /*  default  */ 

if  (argc  >  1)  sscanf  (argv[l]  . '"/.d"  .fcnpts) ; 
if  (argc  >  2)  sscanf  (argv [2]  .'"/.d"  ,fti) ; 
srand(i) ; 

norm  =  1 . / (MAXINT ) ;  /*  machine-dep.  const.,  to  normalize  */ 

for  (i  =  1;  i  <=  npts;  i++)  { 

s[i]  =  norm*rand();  /*  random  by  default  */ 
x[i]  =  norm^randO ; 
y[i]  =  norm^randO ; 

/*  get  numbers  from  stdin  until  blank  line,  then  random  */ 
if  (getmore)  { 

printfC'Map  '/.d.  Enter  scale,  x,  y:  ”,i); 
gets(line) ; 
if  (line[0])  { 

sscanf  (line ,  "V.lf7.1f7.1f"  ,s+i  ,x+i  ,y+i) ; 

if (s[i]<0. )s[i]  =  -s[i];  /*  enforce  0  <=  s  <  1  */ 

while(s[i]>=l . )  s[i]  *=  0.1; 

> 

else  getmore  =0; 

> 

printfC'Map  */,d:  s  =  '/,f  ,  x  =  */,f ,  y  =  */,f\n" ,  i,s[i]  ,x[i]  ,y  [i] ) ; 
> 

return(npts) ; 

> 
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FIGURE  CAPTIONS 

Figure  1.  An  attractor  of  an  IFS  is  shown  with  its  envelope  of  three  disks,  as 
computed  by  the  direct  algorithm.  (This  IFS  uses  affine  maps,  with  sj  =  =. 
$2  =  5.  «3  =  5.  Xj  =  (0,0).  x-2  =  (4.0).  and  x3  =  (0.3). 

Figure  2.  For  an  equilateral  Sierpinski's  Triangle  (where  tr,(x)  =  +  -;x, !. 

the  r,  =  d,j  by  this  method:  in  this  particular  case  the  radii  could  be  half  as 
large. 

Figure  3.  Same  s,  and  x,  as  the  previous  figure,  but  here  tr,(x)  =  — ^x-e  -,x, 
ill  this  case  the  r,  found  above  are  minimal. 

Figure  4.  A  one-dimensional  attractor  constructed  from  the  following  or¬ 
dered  data:  =  3.  *2  =  5.  *3  =  =•  d\2  =  1.  and  <713  =  *  The  r,  are  minimal 

for  such  attractors. 

Figure  5.  Barnsley's  fern  [4]  and  its  envelope  (see  text). 
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